arXiv: 1508.0057lvl [cond-mat.mes-hall] 3 Aug 2015 
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It is well-known that a non-vanishing Hall conductivity requires time-reversal symmetry breaking. 

However, in this work, we demonstrate that a Hall-like transverse current can occur in second-order 
response to an external electric field in a wide class of time-reversal invariant and inversion breaking 
materials, at both zero and twice the optical frequency. This nonlinear Hall effect has a quantum 
origin arising from the dipole moment of the Berry curvature in momentum space, which generates 
a net anomalous velocity when the system is in a current-carrying state. We show that the nonlinear 
Hall coefficient is a rank-two pseudo-tensor, whose form is determined by point group symmetry. We 
discus optimal conditions to observe this effect and propose candidate two- and three-dimensional 
materials, including topological crystalline insulators, transition metal dichalcogenides and Weyl 
semimetals. 
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Introduction — The Hall conductivity of an electron sys¬ 
tem whose Hamiltonian is invariant under the time re¬ 
versal symmetry is forced to vanish. Crystals with suffi¬ 
ciently low symmetry can have resistivity tensors which 
are anisotropic, but, Onsager’s reciprocity relations [1] 
force the conductivity to be a symmetric tensor in the 
presence of time reversal symmetry, implying that, when 
the electric field is along its principal axes, the current 
and the electric field are collinear, at least to the first 
order in electric fields. However, this constraint is only 
in linear response theory and does not necessarily enforce 
the full current to flow collinearly with the local electric 
field. 

In this paper we study a special type of such non-linear 
Hall-like currents. We will demonstrate that metals with¬ 
out inversion symmetry can have a non-linear Hall-like 
current that arises from the Berry curvature in momen¬ 
tum space. The conventional Hall conductivity can be 
viewed as the zero order moment of the Berry curvature 
over occupied states, namely, as an integral of the Berry 
curvature within the metal’s Fermi surface in momen¬ 
tum space. The effect we discuss here is determined by 
a pseudo-tensorial quantity that measures a first order 
moment of the Berry curvature over the occupied states, 
namely its dipolar distribution, and hence we call it the 
Berry curvature dipole. This nonlinear Hall effect has 
a quantum origin arising from the anomalous velocity of 
Bloch electrons generated by the Berry curvature [2] , but 
it is not expected to be quantized. 

In a time reversal invariant system, the Berry curva¬ 
ture is odd in momentum space, fl a {k) = — f2 a (— k), and 
hence its integral weighed by the equilibrium Fermi dis¬ 
tribution is forced to vanish, because states at k and —k 
form Kramers pairs and are equally occupied. However, 
the second order response is determined by the integral 
of the Berry curvature evaluated in the non-equilibrium 
distribution of electrons computed to first order in the 
electric field. Since the non-equilibrium current-carrying 


distribution is not symmetric under k —► — k, the integral 
of the Berry curvature weighed by it can be finite, lead¬ 
ing to a net anomalous velocity and hence a transverse 
current. 

Our study builds upon a seminal work by Moore and 
Orenstein [3], which predicted a DC photocurrent in 
quantum wells without inversion symmetry due to the 
anomalous velocity associated with the Berry phase. The 
quantum nonlinear Hall effect presented here can be re¬ 
garded as a generalization of this effect. We predict that 
an oscillating electric field can generate a transverse cur¬ 
rent at both zero and twice the frequency in two- and 
three-dimensional materials with a large class of crystal 
point group symmetries. In particular, the second har¬ 
monic generation is a distinctive signature that may fa¬ 
cilitate the experimental detection of the quantum non¬ 
linear Hall effect. Additionally, the effect does remain 
finite in the dc limit of the applied electric field. 

General theory — The electric current density is given 
by the integral of the physical velocity of the electrons, 
v a , weighed by their occupation function, f(k): 

ja = -e [ f(k ) v a . (1) 

Jk 

For simplicity we imagine a single band system but allow 
it to be two- or three-dimensional: f k = f d d k/(2n) d . 
The velocity contains two contributions, namely, the 
group velocity of the electron wave and the anomalous 
velocity arising from the Berry curvature [2] (h = 1): 

Va — b) a e{k J T saba^lbkc: (2) 

where e and fib are the energy dispersion and the Berry 
curvature of the electrons in question: 

fla = ^abc^b-^c: A c = Z (ilk\d c (u/j). (3) 

Within the Boltzmann picture of transport, the canonical 
momentum of electrons changes in time in response to the 
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external electromagnetic fields. In the absence of external 
magnetic fields the change of momentum is: 

k c = —eE c (t). (4) 

where E c (t) = $t{£ c e zult }, with £ c £ C, is the driving 
electric field which oscillates harmonically in time but is 
uniform in space. In the relaxation time approximation 
the Boltzmann equation for the distribution of electrons 
is [4]: 

- erE a d a f + rd t f = fo~ f, (5) 


where fo is the equilibrium distribution in the absence 
of external fields. We are interested in computing the 
response to second order in the electric field, hence we 
expand the distribution up to second order: / = 5ft{/o + 
/i + / 2 }, where the term /„ is understood to vanish as 
£ n . One finds a recursive structure: 


A ruj jiujt fuj 

= J i e > Ji 


eT£qd a f 0 
1 + iiOT ’ 


h 


r 0 | n2uj 2lUJt rO 
J2 ' J2 ° i J2 


(er) 2 £* a £ b d ab f 0 
2(1 + iujt) ’ 


f2 uj 

J 2 


(er ) 2 £ a £ b d ab fo 
2(1 + *wr)( 1 + 2 iwr) 


( 6 ) 


By writing the current as j a = + f u e 2lult } 1 one 

obtains: 


fa = \ j £abcttb£*fi - e J^f$d a e(k), 
ff = j J e ab A£ c f? -ef ftdae(k). 


The term j® describes a rectified current while the 
term j 2u describes a second harmonic generated by the 
non-linearity. The second terms that appear in both 
expressions above are completely semiclassical and do 
not require the presence of Berry curvature. However, 
within the approximation of a constant scattering time, 
one finds that these non-linear terms are proportional to 
the integral of a three-index tensor, d a e(k)d bc fo(k), which 
is odd under time reversal and hence they are forced to 
vanish. Therefore the only surviving terms are those as¬ 
sociated with the Berry curvature. 

By writing j° = Xabc£b£* ■ = Xab c £b £ c , one has 

that: 

e 3 T f 

Xabc ^ a dc 0 , : C / {^bfo^^d- (8) 

2(1 +IUJT) J k 

The presence of the factor d b fo will gurantee that only 
states close to the Fermi surface will contribute to the 
integral in the low temperature limit, so that this re¬ 
sponse is a Fermi liquid property [7]. Equation (8) can 
be rewritten as follows: 

e 3 r r 

Xabc = ZadcdXF j : \ / fo {^bkld)- (9) 

2(1 + Vjjt) J k 


This expression (9) for the nonlinear conductivity ten¬ 
sor, Xabc j is the first main result of this work. It shows 
that Xabc is proportional to the dipole moment of the 
Berry curvature over the occupied states, defined as: 


Dab = [ fo (d a S2 b ). (10) 

Jk 

It is interesting to note that this tensor is dimension¬ 
less in three dimensions. At frequencies above the width 
of the Drude peak lot 1 and below the interband 
transition threshold, the prefactor in Xabc becomes in¬ 
dependent of the scattering time, so that Xabc is a direct 
measure of quantum geometry of the Bloch states in mo¬ 
mentum space. In the dc limit or for linearly polarized 
electric fields, the Berry curvature dipole term always 
produces a current that is orthogonal to the electric field 
j a E a = 0 [24], 

To close this section, we wish to remark that there 
exist additional second order corrections to the current 
that arise from modifications to Eq. (2) that are intrinsic 
to the band structure, containing no powers of the scat¬ 
tering time r [5], however, these contributions vanish for 
time reversal invariant systems. Other type of rectifica¬ 
tions might arise in systems with an inversion asymmet¬ 
ric scattering rate, namely one for which the scattering 
from k to k' has a different rate than that from —k to 
—k', which produces a kind of ratchet effect [6]. These 
semiclassical Berry-phase independent contributions are 
distinguished from the quantum nonlinear Hall effect dis¬ 
cussed in this work because they are expected to scale as 

T 2 . 

Berry curvature dipole in three dimensions — Let us ex¬ 
plore the constraints imposed by crystal point symme¬ 
tries on the Berry curvature dipole moment tensor D ab . 
A point symmetry is described by an orthogonal ma¬ 
trix S. Because the Berry curvature is a pseudovector, 
the Berry curvature dipole transforms as a pseudoten¬ 
sor. Hence, crystal symmetries impose constraints of the 
form: 


D = det (S)SDS t . (11) 

To determine which components of this tensor are non¬ 
zero it is convenient to decompose it into symmetric and 
antisymmetric parts: D ± = (D ± D T )/ 2, which trans¬ 
form independently under symmetry operations. The 
antisymmetric part of a pseudotensor transforms as a 
vector, as can be verfied from Eq. (11). The components 
of this vector can be taken to be d a = e a b c D bc / 2. There¬ 
fore for it to be non-zero the crystal must have a polar 
axis. From the 32 crystallographic point groups, 10 allow 
for a polar axis, namely {C n ,C nv } with n = 1,2, 3,4,6. 
The vector d a will be oriented along such axis. The con¬ 
tribution to the current from this antisymmetric part can 
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be written in vector notation as: 

-£* x (d x £), 


r = 


e 3 r 


2(1 + iwr) 


= 9n ,. / x (d x £). (12) 

2(1 + jwt) 


Let us now determine which crystals allow for a non¬ 
zero symmetric part D + . We require the crystal to be 
inversion asymmetric for otherwise the Berry curvature 
would be identically zero due to the time reversal symme¬ 
try. Any real symmetric matrix can be diagonalized and 
has a real spectrum. Let us denote its eigenvalues and 
eigenvectors by Si, e, respectively: D + = X^=i • 

All inversion asymmetric crystals without left-handed 
symmetries allow for D + to be non-zero, but might gen¬ 
erally impose constraints on its eigenvectors to lie along 
the principal symmetry axis and some of its eigenval¬ 
ues to be degenerate, much in the same way they con¬ 
strain an ordinary tensor. Such non-centrosymmetric 
crystal point groups without left-handed symmetries are 
{O, T, C x ,C n , D n } with n = 2,3,4, 6 . 

However under left-handed symmetries (detS = — 1) 
the transformations of D + differ from those of an ordi¬ 
nary tensor. Equation (11) implies that under a left- 
handed symmetry operation the spectrum goes to minus 
itself: {£i, 82 , <$ 3 } —» {—5i, — 82 , — <$ 3 }- Therefore the only 
way that it remains invariant as a set, is if it has the form: 
{<5i, 62 , £ 3 } = {<5, 0,— 5}. In such case the eigenvectors 
would be forced to transform as Se x = ±e 3 , Se^ = ±ei, 
Se 2 = ±e 2 - Therefore, any crystal with a left-handed 
symmetry and an n-fold rotation axis with n > 3 will 
force the tensor D + to identically vanish, since such n- 
fold rotation would additionally force the eigenvectors 
contained within the invariant plane to be degenerate. 
For a mirror symmetry, the null eigenvector has to be 
parallel to the mirror plane, and the eigenvectors with 
opposite eigenvalues must be at 7 r /4 angles from the mir¬ 
ror plane, so that they are swapped under the mirror op¬ 
eration. Therefore, the only non-centrosymmetric crys¬ 
tals with mirror symmetries that allow for a non-zero D + 
can be seen to be C Xv and Civ [25]. For C Xv symmetry 
D + has two independent parameters which can be taken 
to be the positive eigenvalue and the orientation of the 
null eigen-vector within the mirror plane. For Ci v there 
is only one independent parameter, which can be taken 
to be the positive eigenvalue, since the null eigenvector is 
forced to lie along the rotation axis. In addition the crys¬ 
tal point group S 4 , which contains a single left-handed 
four-fold roto-reflection symmetry, allows for a non-zero 
D + , whose null eigenvector is forced to lie along the ro- 
toreflection axis. D + has two independet parameters for 
S 4 , which can be taken to be the positive eigenvalue and 
the orientation of the corresponding eigenvector within 
the roto-reflection plane. 

Berry curvature dipole in two dimensions —In a two- 
dimensional crystal the Berry curvature behaves as 


a pseudoscalar (only the component perpendicular to 
the plane is non-zero), hence the Berry curvature 
dipole behaves as a pseudo-vector contained in the two- 
dimensional plane: 


D a = f fo (3„n,). (13) 

Jk 


This vector has units of length. Therefore symmetry con¬ 
straints are more severe in two-dimensions. In fact, the 
largest symmetry of a 2D crystal that allows for a non¬ 
vanishing Berry curvature dipole is a single mirror line 
(i.e. a mirror plane that is orthogonal to the 2D crystal), 
which would force D a to be orthogonal to it. In vector 
notation the current can be written as: 


r 




2(1 + iuir) 


2(1 + iuir) 


zx£*(D-£), 
z x£(D-£). 


(14) 


The presence of a single mirror symmetry would force 
the linear conductivity tensor to have its principal axes 
aligned with the mirror line. Consequently, according to 
Eq.(14), when the driving electric field is aligned with 
the direction of the Berry curvature dipole vector, D , all 
the current that flows orthogonal to it would arise solely 
from the Berry curvature dipole term. 

Candidate materials — Berry curvature is often concen¬ 
trated in small regions in momentum space where two 
or more bands cross or nearly cross. Therefore, Dirac 
and Weyl materials are excellent candidates to observe 
the quantum nonlinear Hall effect predicted in this work. 
Moreover, since this effect requires the dipole moment 
of Berry curvature, it is advantageous to choose low- 
symmetry crystals with tilted Dirac or Weyl point (see 
below). We propose three classes of candidate materials: 
topological crystalline insulators, two-dimensional tran¬ 
sition metal dichalcogenides, and three-dimensional Weyl 
semimetals. 

The surface of topological crystalline insulator (TCI) 
materials hosts massless Dirac fermions protected by mir¬ 
ror symmetries [ 8 , 9]. In particular, the [001] surface 
of TCIs SnTe, Pbi_ x Sn,„Te and Pbi-^Sn^Se, hosts four 
massless Dirac fermions [10] protected by two mirror sym¬ 
metries. Pairs of Dirac cones with spin-momentum lock¬ 
ing are located near the X points of the surface Brillouin 
zone, forming a Kramers pair. At low temperatures the 
surface undergoes a structural transition into a ferroelec¬ 
tric state and one of the mirror symmetries is sponta¬ 
neously broken [ 11 , 12 ], while the other remains intact. 
As a result two of the surface Dirac cones become mas¬ 
sive, while the other two remain massless [13] (see inset 
in Fig. 1). Since the remaining massless Dirac points 
have vanishing Berry curvature, it is sufficient to con¬ 
sider only the contribution to the Hall current from the 
two Dirac points that become massive in the distorted 
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crystal structure. They acquire Berry curvatures of op¬ 
posite signs, because they are mapped into one another 
by time reversal symmetry. The low energy Hamiltonian 
for the massive Dirac point located at momenta ±A away 
from X\ is given by: 

H s a = v x k x a y - sv y k y a x + sak y + /3a z . (15) 

where s = ±1. /3 is the size of the gap opened by the 
ferroelectric distortion. This low energy theory coincides 
with that previously considered in the literature [10, 11, 
13], except for the term proportional to a which produces 
a finite tilt in the Dirac cones. This tilting effect, which is 
allowed by symmetry and has been observed in ARPES 
studies [14], is required for a nonzero Berry curvature 
dipole (see below). 

The dispersion relation for the titled massive Dirac 
electron is found to be: 

£ s (k) = saky + sign(/z)(/3 2 + v 2 x k 2 x + v 2 /c. 2 ) 1/2 , (16) 

where /i > 0 (yi < 0) for conduction (valence) band. The 
Berry curvature can be found, from Eq. (3), to be: 

= sign(^) sv x v y (3 

s 2 (/3 2 + v 2 x k 2 + i;2fc2)3/2' 1 j 

At zero temperature the Berry curvature dipole is: 

D a = [ d b n. (is) 

This integral can be computed by performing first an 
area-preserving-transformation: k' y = \Jv y /v x k y , k! x = 
\/vx/vyk X i and by noting that the Fermi surface is an 
ellipse in the primed coordinates, namely, that the equa¬ 
tion e s (k ) = yi implies that: k x / q 2 + ( k' y + sko) 2 / 7 ^ = 1 . 
Where -j x = "f/v, j y = 'y/y/v 2 — a' 2 , kg = yia 1 /{v 2 — 
a' 2 ), v = y/Va^Jy, 7 = vV 2 + {v 2 - a , 2 )fc 2 _ ^ a > = 
a^/v x /v y . The condition v 2 > a! 2 is equivalent to 
v 2 > a 2 and is needed for the stability of the Dirac 
cones in the presence of the tilt term. The condition 
yi 2 + ( v 2 — a , 2 )A:o > /3 2 states that the chemical poten¬ 
tial is outside the gap, so that there is a finite density of 
massive Dirac fermions. 

The surviving mirror symmetry, that takes k y —> —k y , 
dictates that only the y-component of the Berry curva¬ 
ture dipole is non-zero, and is found to be: 

n _ 3u 2 n/3a|p|(l + u 2 ) /in ^ 

u y — r c./n i 

[^ 2 (l + u 2 )(l + 2 u 2 )-u 2 p 2 ] ' 

where u = a'/y/v 2 — a ’ 2 and n = /, | d 2 k/( 2 n ) 2 = 

Ixly/^ is the absolute value of the carrier density in 
each of the massive Dirac cones [26] . Each of the massive 
Dirac cones has identical contributions to D y , giving rise 
to a factor of 2 already included in Eq. (19). This dipole 



FIG. 1: (color online) Berry curvature dipole dependence on 
chemical potential /r. The upper inset is the surface Brillouin 
zone of the TCI SnTe or (Pb,Sn)Se. The blue arrow indicates 
the direction of the ferro-electric distortion. The lower inset 
illustrates the Brillouin zone of monolayer TMDC with Dirac 
points shifted away from K and K' by shear strain along 
the directions indicated by the blue arrows. Red circles and 
dashed lines indicate massive Dirac points and the surviving 
mirror symmetry respectively in both insets. 

is orthogonal to the ferroelectric displacement direction 
in the convention we have chosen. The Berry curvature 
dipole has the same sign for electrons and holes in this 
system and vanishes when the chemical potential is in 
the gap of the massive Dirac fermions. The typical scale 
of D y for SnTe TCI is Ha/j3 ~ 3nm, where we used a 
Fermi velocity of v x « v y «4x 10 5 m/s [23], /3 ~ lOmeV 
and a = 0.1v x . The behavior of D y is depicted in Fig. 1. 

Another candidate 2D materials to observe the quan¬ 
tum non-linear Hall effect are monolayer transition-metal 
dichalcogenides (TMDC). Their large spin-orbit-coupling 
and lack of an inversion center produces substantial local 
Berry curvatures [15, 16]. The C 3v symmetry of these 
crystals would force the Berry curvature dipole to van¬ 
ish. However, uniaxial strain can reduce this symmetry 
so that a single mirror operation survives, in which case 
the effect is allowed. In fact, two copies of the model of 
Eq. (15), each with a different gap, can describe the states 
near charge neutrality within a k-p model [15], and when 
the shear strain is applied along high-symmetry lines (see 
inset of Fig. 1). s — ± would label valleys K and K’ in 
this case. The anisotropic velocity term parametrized by 
a would be proportional to the shear strain, much in the 
same way as in strained graphene [17]. For TMDCs one 
obtains a scale Ha/fi ~ 0.2A, using a Fermi velocity of 
v « 4.5 x 10 5 m/s, a gap j3 « 1.5eV and a = O.lu. 

Last but not least, the Berry curvature dipole induced 
non-linear Hall effect should be present in a very large 
class of three-dimensional non-centrosymetric crystals. 
An interesting candidate are the recently discovered Weyl 
semimetals in the TaAs material class [18-21]. The crys- 
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tal structure of these materials is non-centrosymmtric 
and has a polar axis, which allows the quantum non¬ 
linear Hall effect as described by Eq. (12). When tilted, 
a Weyl point generates a singular configuration of Berry 
curvature, with a finite dipole moment whose magnitude 
can be easily estimated from band structure calculations. 
In addition, other polar materials such as BiTel with 
a strong Rashba-type spin-orbit coupling [22] may also 
have large Berry curvature dipole moments. These three- 
dimensional Weyl and Rashba materials provide promis¬ 
ing platforms for the observation of the quantum nonlin¬ 
ear Hall effect. 
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